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We propose that the superconducting pairing symmetry 
of organic superconductors k-(BEDT-TTF)2X can be deter- 
mined by measuring the position in momentum space of the 
incommensurate peaks of the spin susceptibility. We show 
that the positions of these peaks are uniquely determined by 
the symmetry of the superconducting state and the geometry 
of Fermi surface. Based on the effective Hubbard model on 
an anisotropic triangular lattice, we demonstrate the differ- 
ent incommensurate patters of spin responses for d x 2_ y 2- and 
d xy -\ike pairing states. In addition, we predict the nonexis- 
tence of the spin resonance mode in this compound, which is 
an important feature in high-T c superconductors. 

PACS number: 74.70.Kn,74.20.Rp 

Layered organic superconductors based on the (BEDT- 
TTF) molecule are rich examples of strongly correlated 
electron systems in low dimensions . Among them, k- 
(BEDT-TTF) 2 X(k-(ET)) has attracted considerable at- 
tention, because they share many similarities to cuprate 
superconductors , such as the proximity of its super- 
conducting(SC) phase with the antiferromagnetic(AF) 
insulating phase which suggests that it may be an uncon- 
ventional superconductor. Understanding its SC mecha- 
nism could therefore provide important insight into the 
origin of the unconventional superconductivity. As an 
essential step for this process, intensive experimental 
and theoretical investigations have been carried out to 
identify its SC pairing symmetry, but no consensus has 
been achieved yet ||. The NMR Q| and millimeter- 
wave transmission experiment || suggests d x 2_ y 2 sym- 
metry (A different interpretation has been proposed for 
the latter result ||). However, the recent thermal con- 
ductivity and scan tunneling spectroscopy measure- 
ments ^| predict that it is more like d xy than d x 2_ y 2. 
Moreover, a nodeless s-wave SC state is also suggested 
in recent specific- heat measurements ||. Theoretically, 
the d x 2 _ y 2-\\ke pairing symmetry is favored in the model 
in which the superconductivity is mediated by strong AF 
fluctuations (H]|n]] and also in a model of small-g phonon 
mediated superconductivity Jl2fl , However, it is shown 
that the d xy -\ihe superconductivity is also likely and can 
dominate over d x 2_ y 2 wave, when a more realistic model 
is considered [jL3| . 

In this letter, we propose an orientational-dependcnt 
method aiming at a direct determination of the gap 



anisotropy in k-(ET). It is to measure the g-space posi- 
tion of the incommensurate(IC) peaks of the dynamical 
spin susceptibility by inelastic neutron scatterings. Our 
proposal is based on: i) it is believed that AF spin fluc- 
tuations are strong in this material Jl],|2|,|| ; h) the central 
issue at problem about the gap symmetry is what the 
direction of the gap nodes is or if there are gap nodes. 
The underlying physics is that the particle-hole excita- 
tions, which contribute to the spin susceptibility, are only 
available for the node-to-node transitions in the SC state, 
if the excitation energy is lower than 2Ao(Ao the magni- 
tude of the SC gap). Therefore, the g-space position of 
spin responses is uniquely determined by the node posi- 
tion and the geometry of the normal state Fermi surface. 
The Fermi surface in k-(ET) has been determined by var- 
ious experiments ||] , so the node position can be reliably 
determined by measuring the g-space structure of spin 
excitations. This method may be useful specially for k- 
(ET), because the nodes of the d x 2_ y i-wave pairing are 
different from those of d xy -wa,ve pairing by rotating the 
former 7r/4 along the k z axis. 

The above argument is illustrated by a numerical cal- 
culation based on the Hubbard model on an anisotropic 
triangular lattice at half filling, which is adopted as the 
simplest model for k-(ET) pjp^. We find that the spin 
response is incommensurate in the normal state due to 
the nesting of Fermi surface. In the SC state, the IC 
peak appears at (0.97T, 0.97r) for d x 2_ y 2 wave, but around 
(0.357T, 0.357r) for d xy wave. For an isotropic s-wave, the 
peak is suppressed highly and its structure is quite sim- 
ilar to that in the normal state. We also investigate 
the frequency dependence of dynamical susceptibility at 
the AF wavevector Q = (tt,tt) and no resonance peak 
is found in the reasonable range of parameters. It con- 
trases with the case in high-T c cuprates, where a reso- 
nance peak is observed at Q for YBaCuO, BiSrCaCuO 
and TIBaCuO |T§. 

The model Hamiltonian is p4| , 

H = -t 4, - c j>-* / H c l<r c 3,<r + U Yl ni >1 ni ^ 

<ij>,cr <ij>' ■<? i 

(1) 

where the notation is standard except that < ij >' de- 
notes the next-nearest-neighbour bond only in one di- 
rection. Its dispersion is = — 2t(cosk x + cosk y ) — 
2t' cos(k x + k y ) — fi. According to the band structure 
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calculations @, we take if /t « 0.7 and U » 7t. The 
chemical potential /i = 0.7f is determined by fixing the 
electron density to be at half filling. The Fermi surface 
is shown in Fig.l. 

In the SC state, the bare spin susceptibility Xo(q> w) = 
Xo(q, iu n - > uj + iS) is given by, 

Xo(q, iw n ) = -T^ G(k, icj m )G(k + q, iu>„ +iui m ) 

U..TTI 

+ F(k,iuj m )F{k + (i,iuj n + iuj m ), (2) 

where G and F are the normal and anomalous Green's 
functions, 

_ + e k A k 

-e»- A*' (^„)2_ e 2_ A 2' W 

with A,t the SC pairing order parameter. According to 
the recent tunneling experiment Q, the maximum gap 
value A is taken to be 3.0meV=0.04t. 

Bulut et ai, |p"7| have shown that a random-phase 
(RPA) approximation for the spin susceptibility with a 
renormalized (reduced) Coulomb coupling U — 2t gives a 
fit to the Monte Carlo data for a 2D Hubbard model of 
U = At. Thus, we consider the renormalized spin suscep- 
tibility in the RPA approximation which is given by, 

x(q>u) = xo(?,^)(i - Uxoiq^Jr 1 (4) 

where the renormalized vertex U instead of the bare U 
has been used. The system will be instability to an AF 
Neel state if the bare U is used, which is not consis- 
tent with the metallic normal state considered here. This 
is due to the artifact of the RPA approximation, which 
overestimate the AF correlation as was found in the cases 
of the cuprate system |T^] and also organic materials [jl9| . 
Considering the conclusion of Bulut et al. [jl7| and that 
the instability to an AF Neel state occurs when U > 3.4i 
, we choose U = 3t here. 

Four different symmetries of the pairing state are in- 
vestigated here: i) = Ao(cosfe x — cosfc a )/2(pure 
d x 2_ y 2 wave); ii) Afe = Ao sin k x sin k y (pure d xy wave); 
iii) Afc = 0.4A (cos fc-j, — cosk y ) + 0.2A (mixed d x 2_ y 2 
wave); and iv) Afe = 0.8A sin k x sink y ± 0.2 A (mixed 
d xy wave). We note that one should differentiate +0.2A 
and — 0.2A in case iv) due to the anisotropy of the Fermi 
surface. This will be discussed below. 

Firstly, the momentum(q) dependences of the spin sus- 
ceptibility in the normal state at T = O.OIi are presented 
in Fig. 2(a) and (b) for id = 0.02t and ui — 0.2t, respec- 
tively. An obvious and also common feature for both 
cases is that the spin response is incommensurate. The 
IC peaks form two parallel walls near Q = {%, w). They 
come from the excitations with the transition wavevector 
1 connecting the nesting pieces of the Fermi surface as 
shown in Fig.l. 



The spin response in the SC state with d x 2_ y 2 sym- 
metry for u> = Q.02t is shown in Fig. 3(a). The quasi- 
particle spectrum now has gap except the four nodes 
(±fco,±fco) a l° n g the diagonal direction of the Brillouin 
zone. For to = 0.02t, which is lower than the maximum 
SC gap A = 0.04i, the excitations away from the nodes 
are suppressed and only those from nodes to nodes sur- 
vive. As a result, one can only observe two sharp peak 
at (n ± 0.17T, 7r ± 0.l7r)(The (— k ,k ) to (fco, — fco) exci- 
tations do not give rise to strong peaks due to the lack 
of the nesting of Fermi surface along that direction). We 
have also calculated the result for to = 0.2t and find that 
it is nearly the same as that in the normal state, because 
the quasiparticle excitations are nearly not affected by 
the opening of the SC gap. 

We show the results with d^-wave symmetry in 
Fig. 3(b) for u) = 0.02t. Compared with the result of 
d x 2_ y 2 wave, an obvious difference is that a new peak 
appears near (0.357T, 0.357r). Moreover, its intensity is 
stronger than the original one near (tt,tv). The origin of 
this peak can be seen from the Fermi surface shown in 
Fig.l. For dzy-wave symmetry, the four nodes are now 
at (±pi, 0), (0, ±pi), where pi is the cross point of the 
Fermi surface with the k x or k y axis. So, the dominat- 
ing particle-hole excitations at low frequencies are the 
excitations corresponding to the transition 2 as shown 
in Fig.l. The reason that this peak is along the diago- 
nal direction instead of the k x or k y direction is that the 
transition wavevector connecting the nesting piece of the 
Fermi surface is along that direction. Now, let us discuss 
the conditions which this new peak can be observed in 
experiments, i.e., which its intensity is stronger than the 
IC peaks around (tt, tt). Firstly, the excitation energy 
for the spin response should be two times smaller than 
the magnitude of the SC gap Ao. Otherwise, the excita- 
tions (transition 1) across the SC gap will overwhelm the 
node-to-node excitations (transition 2) due to the much 
stronger nesting effect of Fermi surface near (tt/2,tt/2). 
We show the imaginary part of the dynamical susceptibil- 
ity along the diagonal direction for different frequencies 
lu = 0.02t, 0.08i and 0.2t in Fig.3(c). When lo is equal to 
2Ao = 0.08£, the intensity of the peak at (0.357T, 0.357r) is 
nearly the same as that near (tt, tt). Above 0.08£, the for- 
mer decreases and the latter increases, but the opposite 
case occurs for u> < 0.08t. Secondly, the damping rate of 
quasiparticles should be lower than Ao due to the same 
reason. Numerical result shown in Fig. 3(d) supports this 
analysis. As for the case of nodeless s-wave symmetry as 
suggested by recent experiments H , we find that the spin 
response is suppressed to be only 5% of that for d x 2_ y 2 
wave and u — 0.02t, due to the fully opened gap, and 
its structure is quite similar to that in the normal state 
shown in Fig. 2. Therefore, the difference in the positions 
of IC peaks can be taken to differentiate the symmetry 
of the pairing states thus far proposed for k-(ET). 

Due to the underlying anisotropic triangular lattice in 
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k-(BT) [jl|, its SC gap symmetry is not a pure d wave. 
We mimic the possible variation of the SC gap nodes by 
the mixed d and isotropic s wave. This is reasonable, 
because the essential point we suggested in this paper is 
the position of the nodes, instead of the gap shape. In 
the case of Afe = cA (cosk x — cosfc y )/2 + (1 — c)A , the 
nodes at (fco,feo) and (— fco, — kg) shift clockwise. There- 
fore the IC peak which is along the diagonal direction for 
pure d x 2_ y 2 wave also shifts clockwise. This can be seen 
clearly from Fig. 4(a) with c = 0.8. Because the shape of 
Fermi surface is symmetric with the axis along the diag- 
onal direction, the ±0.2Ao terms will have a symmetry 
effect. However, the situation is different for a mixed d xy 
wave with a small s-wave component, due to the lack of 
the symmetry of Fermi surface with the axis k x or k y 
which is the node direction of a pure d xy wave. For a 
small positive s-wave component, the node at (pi, 0) will 
shift downwards from the k x axis to the position being 
more near the nesting wavevector (0.357T, 0.357r). There- 
fore, the IC peak at (0.357T, 0.357r) is enhanced and has 
an much more overwhelming intensity than those peaks 
around {%, n), as shown in Fig. 4(b). But for a negative s- 
wave component, the node at (pi, 0) will shift oppositely, 
i.e., move away from the nesting wavevector, and the 
peak at (0.357T, 0.357r) is strongly suppressed(Fig.4(c)). 
However, this mixed gap symmetry can still be distin- 
guished from the pure or mixed d x 2_ y 2 wave. Because 
four peaks around (7r,7r), instead of two peaks as in the 
pure or mixing d x i_ y i symmetry, should be observed as a 
result of the suppression of the peaks along the diagonal 
direction where the gap is fully opened now. 

Finally, we discuss the w-dependence of the spin re- 
sponse, which is shown in Fig. 5 for (a) d x 2_ y 2 wave and 
(b) d xy wave, respectively. One can see that no reso- 
nance peak is found at momentum q = (-7T, 7r)(also for q = 
(0.357T, 0.357r) in the case of d xy wave). In the particle- 
hole channel, a spin resonance is expected if we have 1 — 
URexoiQ, w) =0 and negligibly small Irnxoil, w) fhSl^of . 
This satisfies in high-T c cuprates [^8|,^0), due to a sharp 
step-like rise of the bare Im%o &t its gap edge, which 
comes from the extended van Hove singularity near (0, ir), 
then a logarithmic singularity in Rexo via the Kramers- 
Kroenig relation. However, it does not satisfy here as 
see from the inset of both figures 5(a) and (b), where 
the largest value of Rexo is about 0.249 and 0.262, re- 
spectively. The reason is that the step-like rise is very 
small or lack in Im^Oj because of the different energy 
band structure in k-(BT). So, no resonance peak can be 
seen even for the largest value U — 3.4£ as shown in 
Fig. 5 (a). The only condition for the occurence of a reso- 
nance peak here is limited to be U — 4.0i±0.1i for d x 2_ y 2 
wave at (tt, 7r), but it is beyond the reasonable parameter 
range as discussed above. The above result is in sharp 
contrast with the w-dependence of the spin responses in 
high-T c cuprates and should have important insight on 
the differences of the superconducting properties in both 



materials pq . 

In summary, we propose an orientational-dependent 
method to determine the SC pairing symmetry in n- 
(BEDT-TTF) 2 X based on the g-space position of the 
incommensurate peaks of the spin susceptibility. It is 
demostrated that the incommensurate peak positions at 
low energies are uniquely determined by the pairing sym- 
metry and the geometry of Fermi surface. Different pat- 
terns of the spin responses in the momentum space are 
shown for d x 2_ y 2 wave, d xy wave, mixed d x 2_ y 2 wave and 
mixed d xy wave. In addition, we predict the nonexistence 
of the spin resonance mode in this material. 
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FIGURE CAPTIONS 



Fig.l Fermi surface in the first and the extended Bril- 
louin zone. The arrows show particle-hole excitations 
with nesting wavevectors (0.97T, 0.97r) (labeled as 1) and 
(0.357T, 0.357r) (2), respectively. 

Fig. 2 Momemtum dependence of lmx(q, lu) in the nor- 
mal state (temperature T = O.Oli) for (a) lu = 0.02i and 
(b) lu = 0.2t. 

Fig. 3 Momemtum dependence of Im%(q, lu) in the su- 
perconducting state (temperature T = 0.00ft) with lu = 
0.02t for (a) d x 2_ y 2 and (b) d xy . (c) and (d) are the re- 
sults along the diagonal direction with d xy for different 
excitation energies and different quasiparticle damping 
rates, respectively. 

Fig. 4 Momemtum dependence of Imx(q, lu) in the su- 
perconducting state (temperature T — 0.00ft) with lu = 
0.02t for (a) 0.8d x 2_ y 2 +0.2 isotropic s wave, (b) d xy + 0.2 
isotropic s wave, and (c) d xy — 0.2 isotropic s wave. 

Fig. 5 Frequency dependence of lmx(q, lu) for (a) 
d x 2_ y 2 and (b) d xy waves. The insets show the frequency 
dependence of the bare spin susceptibility lmxo(q, lu). 
Solid line is their real part and dashed line their imagi- 
nary part. We have added the imaginary part a constant 
0.23 and 0.22 for d x 2_ y 2 and d xy , respectively, in order 
to plot it in the same figure with their real part. 
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